Self-sustained oscillations in a Large Magneto-Optical Trap 
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We have observed self-sustained radial oscillations in a large magneto-optical trap (MOT), con- 
taining up to 10 10 Rb 85 atoms. This instability is due to the competition between the confining 
force of the MOT and the repulsive interaction associated with multiple scattering of light inside 
the cold atomic cloud. A simple analytical model allows us to formulate a criterion for the instabil- 
ity threshold, in fair agreement with our observations. This criterion shows that large numbers of 
trapped atoms N > 10 9 are required to observe this unstable behavior. 
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A large fraction of the stars in the upper Hertzsprung- 
Russell diagram present pulsations based on an inter- 
play between modulated radiation pressure effects, which 
tends to increase the size of the star, and a collapse 
based on gravitational forces . Instabilities also occur 
in other similar systems such as confined plasmas where 
a long range Coulomb interaction has to be countered 
by a confining force to avoid an explosion of the plasma 
[2j . These systems are of fundamental importance for as- 
trophysics and for controlled fusion and have thus been 
extensively studied in the past. However it is either im- 
possible (in the case of stars) or extremely difficult (in 
the case of confined plasmas) to perform experiments 
to study the full dynamics of such systems where col- 
lective effects play a dominant role. On the other side, 
allowing for adequate rescaling, alternative systems can 
present similar dynamics. A variety of interesting collec- 
tive effects have thus been identified in charged colloidal 
systems [3|. Recently ultra-cold plasmas created by ion- 
izing a cloud of laser cooled atoms became subject to in- 
creased attention p| . Beyond the possibility of studying 
analogous effects as in astro- and plasma physics, sys- 
tems with long range interactions are known to lead to 
non-extensive behavior and appropriate scaling laws are 
needed to predict macroscopic properties. Here we show 
that a large cloud of laser cooled atoms is an adequate 
system to study such collective effects. The radiation 
pressure of the multiply scattered photons in such clouds 
can indeed be related to a long range Coulomb type in- 
teraction We thus suggest an analogy between the 
dynamics of a large cloud of cold atoms, astrophysical 
systems and plasma physics. 

The effect of multiple scattering on the dynamics of 
the atoms is well known in the community of laser cool- 
ing of atoms, as multiple scattering has been a major 
limitation to obtain large phase space densities in cold 
atomic traps. Bosc-Einstein condensation (BEC) in di- 
lute atomic vapors has only been achieved after switching 
off all laser fields and using evaporation techniques [y]. 
More recently, multiple scattering of light in cold atoms 
has been used to study coherent light transport in ran- 



dom media pj. This has led to an investigation of yet 
unexplored regimes, namely the limit of very large num- 
ber of cold atoms in the presence of quasi-resonant light. 
Here we do not focus on the properties of the scattered 
light but on the mechanical effects of this light on the 
atoms. We have observed collective instabilities triggered 
by the repulsive interatomic force arising from multiple 
scattering, and identified a supercritical Hopf bifurcation 
separating the standard stable MOT operation from a yet 
undescribed unstable regime. 

In order to estimate the relevance of plasma physics 
considerations to study multiple scattering of light by 
cold atoms it is worth deriving the equivalent of sev- 
eral plasma parameters for our system. The analogy 
with an 1 jr 2 repulsive Coulomb- type force [5j is obtained 
from evaluating the power scattered by one atom (P sca tt) 
and deriving the intensity incident I2 on a second atom 
via I2 cx P sca tt/(47rr 2 ). The resulting radiation pressure 
force scales as 1/r 2 and one can thus define an effective 
charge q which depends on the absorption cross sections 
and laser intensity and is typically q w 10" " 4 e@. A to- 
tal interaction energy qV — ^ R larger than the ki- 
netic energy fc^T of the particles leads to an increased 
diameter L = 2R of the magneto-optical trap when the 
number N of atoms exceeds m 10 5 . Alternatively the 
Debye length A^> = ^e^ksT /nq 2 above which collec- 
tive effects become important is of the order of 
well below the typical size of a large MOT (several mm). 
Also, in our experiments the corresponding plasma fre- 
quency iVD — \/nq 2 /meo is slightly larger (w 200Hz) 
than relaxation rate of the atomic positions (~ 50Hz). 
We thus expect our cloud to behave as a weakly damped 
plasma. Another interesting quantity is the ratio be- 
tween the nearest neighbor Coulomb interaction and the 
kinetic energy T cb = j^/k B T with a ~ n^ 1 / 3 @|. We 
estimate this quantity to be smaller than unity in our 
system, excluding thus any crystallization. An impor- 
tant aspect of these light induced collective interactions 
is that the effective charge q depends on experimental 
control parameters, allowing for a engineering of the ef- 
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FIG. 1: Fluorescence of part of the MOT during a loading 
sequence. Below a critical number of atoms Nth, the size 
of the atom cloud increases without specific dynamical be- 
havior. Above the threshold Nth the cloud switches to an 
unstable mode characterized by periodic oscillations in the 
partial fluorescence signal. Inserts : Fourier transform of sig- 
nal, with (a) a flat noise in the stable regime and (b) distinct 
oscillations in the unstable regime. 

fective charge which can be modified by orders of magni- 
tude. Finally it might be possible to use the high phase 
space densities of a BEC and thus study strongly cou- 
pled plasma in the degenerate regime as expected in 
neutron stars and white dwarfs [lfjj. 

Our cloud of cold atoms is confined in a MOT using 
laser- induced forces We collect Rb S5 atoms from 

a dilute vapor using six large independent laser beams 
(beam waist 4cm, power per beam P = 30mW) thus 
avoiding the intensity imbalance and feedback mecha- 
nism responsible for the instability of ref^l- Under 
standard operating conditions, the trapping lasers are 
detuned from the F = 3 — ► F' = 4 transition of the D2 
line by 6 — — 3T (T/2tt — 6MHz). A magnetic field gra- 
dient (V-B w lOG/cm) is applied to generate a spatially 
dependent Zeeman shift yielding the restoring force of 
the trap. A repumping laser on the F = 2 — > F 1 = 3 of 
the D2 line is used to control the total number of atoms. 
We thus obtain a MOT with u p to N = 10 10 atoms (di- 
ameter L — 5mm, T = 80/i-fT) [13| . The size and shape 
of the cloud is monitored by imaging the MOTs fluores- 
cence on a cooled CCD. The optical thickness b of the 
cloud at the trapping laser frequency is measured by a 
photodiode. To obtain a time-resolved information on 
the local density of the MOT, we also image a portion of 
the cloud on another photodiode. 

FigEJillustrates the onset of spontaneous self-sustained 
oscillation for a sufficiently large number of atoms. We 
switched on the MOT at t = and monitored the time 
evolution of the fluorescence from a portion of the MOT. 



FIG. 2: Fluorescence of the MOT center. In the unstable 
regime, periodic oscillations appear in the absence of external 
modulations. 



This partial fluorescence signal is roughly describing the 
number of atoms in the observed region. Starting from 
N = at t = 0, the trap fills with a time constant 
t = 1.45s determined by the ambient Rb pressure. Be- 
low a critical number of atoms Nth, the size of the atom 
cloud increases with number of trapped atoms [l^ but 
no specific dynamical behavior is observed. Above the 
threshold Nth the cloud switches to an unstable behavior 
characterized by periodic oscillations in the partial fluo- 
rescence signal. These radial oscillations of the cloud are 
self-sustained in the sense that no external modulation 
of any control parameter is present. 

In the inserts of Fig^are shown Fourier transforms of 
the partial fluorescence. Below the instability threshold 
(insert a) a flat noise background is obtained. In con- 
trast, in an unstable MOT, obtained for a larger number 
of atoms, distinct oscillation frequencies (insert b) ap- 
pear, with higher harmonic components indicating the 
non harmonic oscillation of the signal. 

Indeed, the dynamics in the unstable regime can be 
more complex than a harmonic oscillation, as further il- 
lustrated in Fig[21 where we detect the fluorescence from 
the center of the MOT. A high contrast modulation of 
the center fluorescence is observed in this experiment. 
We can speculate that to the fast phase of decrease of the 
signal corresponds a MOT expansion (decreased density 
at the center), whereas we associate the increasing part 
of the fluorescence to a slower compression phase. We ob- 
served that the precise shape of this oscillation depends 
on the laser beam alignment and on the monitored region 
of the MOT. However, the threshold separating the sta- 
ble from the unstable regime was found to be very robust 
with respect to trap parameters. 

Investigating the MOT at the instability threshold by 
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FIG. 3: Phase diagram of self-sustained oscillation : a 
(8, X7B) cut in parameter space shows the separation between 
the stable (for larger detuning 5) and the unstable regime: ex- 
perimental threshold values (squares) and 1-zone model pre- 
diction (solid line). 

varying 2 of the control parameters of the experiment (de- 
tuning, magnetic field gradient), we can map the phase 
diagram shown in Fig|3(full squares). The solid line cor- 
responds to the theoretical prediction presented at the 
end of this paper. As can be seen, the overall behavior 
is unstable when the trapping laser frequency is brought 
within roughly one natural width from resonance. This 
critical detuning depends here rather weakly on the mag- 
netic field gradient. Note however that the measured 
cloud size and number of atoms at threshold do vary quite 
a bit during this experiment (a factor 5 for N and a factor 
2 for L). In addition, we systematically found an optical 
thickness b « 1 at the instability threshold. However, 
this is clearly not a sufficient condition for the onset of 
instabilities, since b = 1 is also observed in the stable 
region of Fig|3] 

We also carefully monitored the fluctuations of the to- 
tal number of atoms when the MOT operates in the un- 
stable regime, as e.g. for the data in FigQl We found 
these to be below 2%, as in the stable regime. This in- 
dicates that, for a given set of MOT parameters, the 
unstable cloud oscillates at fixed N. 

To further characterize the transition to instability, we 
have checked that the amplitude of the oscillation contin- 
uously grows from zero as the control parameter (detun- 
ing or number of atoms in the experiments performed) 
crosses the threshold value. A Fourier analysis of the 
signal across the threshold showed that the instability 
start at a non zero frequency, which is closely related to 
the natural oscillation frequency of the harmonic trap. 
Furthermore, no hysteresis was observed despite explicit 
investigation. All these findings are consistent with a 
supercritical Hopf bifurcation. 



Already in the stable regime, we observed some 
clear indications that strongly increasing the number of 
trapped atoms affects the way the MOT operates. As it 
is well-known, the MOT inflates when atoms are added 
as a consequence of multiple scattering of light ||. In 
addition to the standard L oc we found 

for large number of atoms N > 10 9 a different scaling 
L oc By monitoring the relaxation of the 

MOT after displacing it from its equilibrium position, we 
observed a cross-over from an over-damped behavior at 
small N (typical for usual MOTs) to an under-damped 
behavior at large N . We interpret this finding as a con- 
sequence of the attenuation of the trapping beams inside 
the cloud, which reduces the friction at the center of the 
cloud. This could be envisioned as a precursor to the in- 
stability. Indeed, we found that just below the threshold 
(i.e. in stable operation), the MOT is systematically in 
the underdamped regime. 

To explain the apparition of this new instability, we 
developed a simple model where the screened compres- 
sion force of the MOT is competing against the repul- 
sive interaction due to multiple scattering of light inside 
the cloud. We stress that this instability is thus qualita- 
tively different from that studied in ref 12], where the use 
of retro-reflected beams introduces the feedback neces- 
sary for the instability. The instability process of ref [12, 
which manifests as oscillations of the center-of-mass of 
the MOT, does not involve the long-range interatomic 
interactions which drives the behavior of our large MOT. 

We propose in the following a very simple 1-zone model 
which exhibits an instability threshold. This model 
amounts to an extremely simplified mean field theory, 
based however on microscopic expressions for the light 
forces acting on the atoms. A more refined approach, 
beyond the scope of this paper, could e.g. involve hydro- 
dynamical approximations [l£|. 

We assume an homogeneous density and the size of the 
cloud L is related to the density n via the total number of 
atoms N: n — N/L 3 . The dynamics along one symmetry 
axis (Ox) of a probe particle located outside of the cloud 
(at position x > R from the trap center, with a velocity 
v) is then governed by the force : 

HkT er h 

F(X,V) - — Sin Ci | 4(j_ Alx _ fc ^)2- 
1 "T p2 

hkv l 

2~ Sinc 1 , 4(S+tJ.x+kv) 2 

J- T p2 

(1) 

This expression relies on the low intensity Doppler 
model for the magneto-optical force (incident on- 
resonance saturation parameter Sj nc ). The first term in 
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this expression is the attenuated force of the laser passed 
through the cloud (with the corresponding Zeeman shift 
\ix and Doppler shift kv), the second term corresponds 
to the non attenuated force of the laser propagating in 
the opposite direction. In absence of the e~ b attenuation, 
these two terms give rise to the standard cooling (via the 
opposite Doppler terms kv) and trapping (via the oppo- 
site Zeeman terms fix) of cold atoms. The last terms is 
the sum of all binary repulsive interactions which, using 
Gauss theorem, yields an 1/r 2 repulsion for a probe parti- 
cle outside the cloud. This term can be understood as the 
radiation pressure originating from the MOT with a to- 
tal radiated power corresponding to attenuation of the 6 
laser beams. Here rj corresponds to the ratio between the 
absorption cross section of the incident laser frequency 
and the inelastically rescattered photons []|. We now 
apply this model at the edge of the cloud (x = R). A 
linear stability analysis, using x = R + Sx e luJt around 
the fixed point F(R, v — 0) = , yields the threshold 
condition for an instability (Jm(cj) = 0) : 

C(6,»,b,R)=e + = (2) 



We find that for our experimental parameters the 
threshold is given with a good approximation by : 

S + fiR w (3) 

It should be stressed that in standard MOTs one usu- 
ally has fiR <C | <5 1 : for \S\ = T and a magnetic field 
gradient VB — 10 G/cm, condition © implies a MOT 
diameter of 8 mm ! It is only with N in the 10 10 range 
that such MOT sizes can be obtained. In this regime, 
the edge of the cloud is now exploring the nonlinear part 
of the magneto-optical force. It can be shown from the 
expression of the force that the threshold condition © 
corresponds to the passage from a positive to a negative 
friction at the edge of the cloud. Thus, a small velocity 
fluctuation is amplified instead of damped and the atoms 
at the edge are kicked away from the center of the cloud. 
The cloud thus expands until its optical thickness drops 
below a certain level, where the MOT is back to its stan- 
dard mode of operation (weak repulsion) and the atoms 
are pushed back toward the center. 

In order to further confront our simple model to the ex- 
periment, we computed the threshold value of the detun- 
ing 5 using the control parameter VB and the measured 
values of the size of the cloud R and its optical thickness 
b at threshold. The result corresponds to the solid line in 
FigED an d gives the correct order of magnitude and be- 
havior for the instability threshold. This gives us further 
confidence about the qualitative validity of our model, as 
the size of the cloud varies by more than a factor 2 along 
the threshold boundary. Furthermore, both the thresh- 
old criterion @ and the simple picture of the unstable 



dynamics near threshold are confirmed by recent results 
from a more involved model|l7|. which includes the inho- 
mogeneous density distribution inside the cloud (N-zone 
model) . 

We have described in this paper the observation of 
self-sustained oscillation in a large cloud of laser cooled 
atoms, arising from the competition between the MOT's 
confining force and the long range multiple scattering re- 
pulsive interaction. This new instability process affects 
the behavior of large MOTs containing more than 10 9 
atoms, a regime increasingly found e.g. in experimen- 
tal setups producing BECs. This observation shows that 
large clouds of cold atoms still present a rich dynamics 
with a variety of yet unexplored regimes. A simple 1- 
zone model has been presented which allows to predict 
the instability threshold and to understand the underly- 
ing physical mechanism. Future possible investigations 
include the forced oscillation regime, the spectroscopy of 
excitation modes in this system, gas-liquid-crystal phase 
transitions in the degenerate regime and feedback mech- 
anisms allowing for stabilization of a large cloud of in- 
teracting particles. Progress on the theoretical aspects 
of the systems described in this paper include exploit- 
ing mean field theory and molecular dynamic simula- 
tions. This should allow for a better understanding of 
the bifurcation observed in our experiment and lead to 
study statistical (thermodynamic) properties across the 
threshold. If the degenerate regime could be reached (or 
inducing similar interactions in a Bosc Einstein Conden- 
sate) a mean field theory based on binary collisions as 
in usual Gross Pitaevskii equation will not be valid due 
to the long range interaction, connecting this system to 
strongly correlated quantum systems. 
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